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Introduction 


Joseph Fourier, whose name forms part of the title of this unit, has already been. 
_ Mentioned in this course, in Unit 12, on heat transfer. The present unit is about a 
different part of his work — one which led to results which have become an 
indispensable part of the theoretical tool kit of applied mathematicians, physicists 
and engineers. 


The fundamental idea is that many functions can be represented by the sum of a 
series of sinusoidal terms. In other words, sinusoidal functions may be considered 
to be the building blocks from which other functions may be built up. The point 
of such a procedure is that in many cases a series of sinusoidal terms is easier to 
deal with than the original function. 


The study of the use of sinusoidal functions to represent other functions is called 
Fourier analysis and is the subject of this unit. 


Study guide 


The sections of this unit are intended to be read in the order in which they appear. 
The first four sections all teach new and assessable material. Section 5 does not 
teach any new material but consists of problems on the material of the rest of the 
unit. You should try as many of these as you have time for without looking at the 
solutions. Then compare your answers with the given solutions, and also look at 
the solutions of any problems you did not attempt— you may well find them 
instructive. 


There is a television programme associated with this unit. You should preferably 
have worked through the first three sections before you watch it but, even if you 
have not done so, it is worth watching the programme rather than missing it 
deliberately. There is no television section in the main text of this unit but there 
are notes about the programme in Appendix 3; you should read these notes before 
watching the programme. There is no tape associated with this unit. 


1 What Fourier analysis is used for 


The Fourier method is to express a function as a sum of sinusoidal terms. This 
representation is particularly appropriate when the function is periodic. In this 
section we shall look briefly at periodic functions and then, through examples, see 
why it is useful to be able to express a general periodic function as a sum of 
sinusoids. | 


1.1 Periodic functions 


A function is said to be periodic if its graph consists of a continually repeated 
pattern, as in Figure 1. We can express this more formally by saying that a fy) 
function f(t) is periodic when 


f(t) =f(t + nT) for all t, forn = +1, +2, +3,.... (1) 


Definition (1) implies that the graph of f(t) is unaffected by being shifted a 
distance nT to the right or left. 


; period 
The positive constant T is called a period of the function. Usually T is taken to be 


the size of the smallest repeating unit in the graph of the function, as in Figure 1, 
but note that we could take any integer multiple of this basic repeating unit as the f(t) 
period since this, too, will satisfy the definition (1). 


When we want to apply the Fourier method to a general periodic function we 
shall need to be familiar with the periodic behaviour of our building blocks, the 
sine and cosine functions. For example, the function sin (2xt/T) has period T 
because 


(b) 
el RSL = sin e = 2nn| = sin la Figure 1. Two periodic 
ie T iy functions 
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so we would expect the function sin (22t/T) to be useful when we approximate 
functions with period T. 


Question 


. 2amt 
Does the function sin , where m is an integer, have T as a period? 


2mm 
Does cos 


Answer 
Yes to both questions, because 


_ 2uam(t+nT) . [2nmt 
sin ————————- = sin + 27mmn 
T 
_ 2nmt ; : 
== Sin 7 because mn is an integer, 
and 
2mmt 


2nm(t + nT) 
ee a =C 


by a similar argument. 


We shall see later how these functions are used in constructing Fourier 
representations of functions with period T. 


1.2. An application of Fourier analysis 


Fourier’s idea of using a sum of sinusoidal terms to represent a periodic function 
turns out to be very useful in many different fields. One of these is the response of 
mechanical systems to any arbitrary periodic input. Since you have some 
knowledge of such systems from Units 7, 8 and 24, the following examples and 
exercises are drawn from that field. 


Example 1 


Assume that the system we want to study can be modelled by the lumped- 
parameter spring-mass system with linear damping shown in Figure 2. (This is 
similar to the systems you modelled in Unit 8.) The force, F, is periodic but non- 
sinusoidal. We want to know how the position, x, of the particle, measured from 
its equilibrium position, will vary with time. How can the Fourier representation 
of F be used to find x? 


Solution 
The equation of motion is 


mX(t) + rx(t) + kx(t) = F(t). (2) 


It could well be very difficult to solve this equation with F in its original form, as 
a general periodic function. But now suppose that we can express F as a sum of 
several sinusoidal terms, say 


F=F, +F,+F3+°°'+ Fy. 


We can use the fact that Equation (2) is a linear second-order differential equation. 
This means that the superposition principle can be applied. Thus if we can find the 
steady-state outputs, x1, X2, ..., X,, that would be produced by the component 
inputs, F',, F2,°-*, F,, acting separately, then we know that the steady-state 
output, x(t), produced by their sum, F, will just be the sum of these: 


x(t) = x,(t) + x2(t) +--+ + x, (t). 


Now, since the terms F,, F2,..., F, are sinusoidal, we have reduced the problem 
of finding x(t) to that of finding the output, x,(¢), produced by a single sinusoidal 
term, F(t). What we need to know is the frequency response of the system, that is 
the way in which the output produced by a sinusoidal input of given amplitude 
depends on the frequency. But we do know how to find the frequency response, 
using the methods of Unit 8. 


In this unit I shall use the results we obtained in Unit 8, but with some changes in 
notation. (The notation I use here is common in the sort of engineering 
applications where Fourier methods are used.) In Unit 8 we represented a general 


Strictly speaking these two 
functions have smallest period 
T/m but, since any integer 
multiple of T/m is also a 
period, the functions also have 
T as a period. 


Figure 2. Model system. Ip is 
the unstretched length of the 
spring. 


Equation (2) can be derived 
by the methods of Unit 8. 


You met the superposition 
principle in Section 2.4 of 
Unit 6. 
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sinusoidal input of angular frequency w by Ag cos (wt + @o). Here I shall use the 
input Fo sin wt, which corresponds to @y = —7/2. I shall write the corresponding 
steady-state output as 


x = Asin(wt + ¢). 


@ is the phase of the output relative to the input, and is often called the phase shift. 


We want to find x for a given input amplitude F, and different values of the input 
angular frequency, w. We can do this by investigating the behaviour of the relative 
response, A/F, and the phase shift, @, for varying values of the angular frequency, 
w. This can be represented by the two graphs sketched in Figure 3, which together 
show the frequency response of the system. The exact shapes of the graphs will 
depend on the amount of damping. Notice that in this example @ is negative for 
all values of w. 


Graphs like those in Figure 3 enable us to read off values of A/Fo and of @ for 
any given angular frequency. So we can find the output, x,, for each sinusoid, F;, 
that makes up the original input. We can then add up all the x; to find the total 
response to the periodic input, F. 


In the above example we used the fact that the system was modelled by a linear 
second-order differential equation, so that we could use the superposition 
principle. More complicated systems may lead to higher-order differential 
equations, but so long as the equations are linear, the superposition principle still 
applies. We define a linear system as a system that can be modelled by a linear 
differential equation. 


We can now generalize the method of Example | to find the response of any linear 
system to a periodic input, F. If we can express F as a sum of sinusoidal terms 
then all we need is a knowledge of the frequency response of the system. This may 
be found. analytically, as in Example 1, or experimentally, by observing the 
response of the system to sinusoidal input forces of various frequencies. (Such 
experiments afe standard practice in many practical situations.) 


Example 2 


A linear system is subject to a periodic input, f(t), which can be expressed as the 
sum of two sinusoids as follows: 

f(t) = 2 sin 5t + 1.5 sin 10t 
where t is time. 


(i) Figure 4 shows graphs of y, = 2 sin 5t and y, = 1.5 sin 10t. Plot f(t) and 
determine its period. 


y 

3 

2 Py a . 

a *, y, =2 sin St 
Fe e | 
fa * ef ad y2 = 1.5 sin 10¢ 
1 me \ * / \ 
ms \ : 7 


Figure 4 


In Unit 8 we used the phase 
lag, Ad, to describe the 
difference between input and 
output. The phase shift is 

= —Adg. 


AlFy 
0 w 
ob 
0 Ww 
— m/2 
— 


Figure 3. Frequency response 
for the model system in Figure 
Pg 


The definition of a linear 
second-order differential 
equation in Unit 6 extends 
naturally to higher orders. 
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(i1) Figure 5 shows the frequency response of a linear system. Derive an 
expression for the response of this system to the input f(t). 


output amplitude 
input amplitude 


1.0 cuceeucceusceccscces sect scsi 


0 : 5 ascgeeee Stascsceance 


Saar 
ie can 


eaagene Ssnee Ress 
HH Sseeeeeeeeeceacceceasaa SeSSG000 G00 CCUS0C0GR00sS5555 Sc00n Hees see 
BESS VZES SAESTESESE ADRES SRSES SERRE REE SS HORSE ENR RE ARREST ERAES SOLERO SHS BROS ERSRSE ROSES EBOBE Dew 
SRBREBR Sa CRs Seaeeeeens seeeeeen Re ees SRS SS SESE SSRSSSSSSS SESS REESE Seees eee t aee8 
SHSBREE HE SRSREPRHES SERERSORAR SE ASPRAR ES GPSAGEREEE PER EG ESR ED SARE CRABS SSERESERSS eS 
SRERRERE V2 CHR AS CRORE SHRRERERAE EBSHARSERE CSRARE DARL DEERE LESH BRIE DAES ES eee h aa 
on Strecrscees 


cans 

{| a ae a8 eaa.—-2 ae rT ESGSsEaSS a t 
mpefetit tt ttt tp CORPS ELE BER RRZERTLOHE EERE ERE AS 
| | @ et ttt ttt SER GRHEREERERRARERRREES 
ae San EEC CEC EEE | BEERS Bun 


Note: the phase shift is 
measured in degrees. 


Figure 5. Frequency response of a linear system. 


Solution 


(i) , 


3 fo Te at Toe 


2 sin St 


3 ok y2 = 1.5 sin 10¢ 
‘ Ve — 2 ; 
: / \ 


\ » Fe 
\ / / 
‘ / f / 


—3 


Figure 6 
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Figure 6 shows both components, y, and y2, of f(t) and their sum, y, which equals 
f(t). The pericd, T, of f(t) is equal to that of the component 2sin 5t; in other 
words 27/5 ~ 1.26. 


(ii) To find the response to f(t) we find the output corresponding to each input 
term separately. 


For the output corresponding to the term 2 sin 5t we need to find the relative 
response and phase shift for the frequency of this input, i.e. w = 5. 


From the top graph in Figure 5 we read the relative response for w = 5 to be 
about 0.93. So the output amplitude for the input 2sin 5t is 0.93 times the input 
amplitude, i.e. 


093 x2 = 19. 


From the lower graph in Figure 5 we read the phase shift, @, corresponding to 
w = 5 to be about — 26°, or —0.45 radians; so the component of the output 
contributed by the input term 2 sin St is approximately 


1.9 sin (St — 0.45). 


(Remember, from Unit 8, that the frequencies of input and output for a linear 
system subject to a sinusoidal input are equal.) 


Similarly, for the output component corresponding to the input term 1.5sin 10¢, 
we read off the relative response and phase shift for w = 10. The relative response 
is about 0.72, so the output amplitude at w = 10 is 0.72 x 1.5 ~ 1.1. The phase 
shift is about —46°, or —0.8 radians. Hence the output component contributed by 
the term 1.5sin 10t is approximately 


1.1 sin (10¢ — 0.8). 


Consequently, by the superposition principle, the total response of the system to 
f(t) is about 


19sin (St — 0.45) + 1.1 sin (10t — 0.8). 
The point that is made in part (i) of the solution is a particular example of a more 
general principle: 


If you take the sum of a number of sinusoidal functions of the form A, sin at, 

A, sin 2mt, A; sin 3at, ---, A, sin kat, so that all the frequencies are integral 
multiples of the lowest frequency, then the sum is periodic and the period is equal 
to that of the component with the lowest frequency (often called the fundamental 

- angular frequency). In other words the function 


f(t) = A, sinwt + Az sin 2mt + A3sin3mt + --:+ A,sinkot 
is periodic and its period is T = 2z/o. | 


Exercise 1 


Suppose the variation with time t (in seconds) of a voltage at a point in an electrical - 
network is given by 


V(t) = 0.01 sin 2t + 0.005sin 4t + 0.002 sin 8¢. 


Determine 
(1) the fundamental angular frequency; 
(ii) the period of the voltage. 


[Solution on p. 23] 


Exercise 2 


Figure 7 shows the frequency response of a linear system. The input of this system is given 
by 


f(t) = sin 4t + 2sin 8t + 3 sin 12t 
where t is the time in seconds. 


(i) What is the period of f(t)? 
(11) Derive an expression for the response, x(t), of the system to the input, f(t). 


[Solution on p. 23 | 
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output amplitude Ft Ee ears ca aia 
input amplitude — FEES Sain ate te Steet 


BERSeeoee — tp eae 
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phase shift, ¢ 
(degrees) 


oe: stitsess 
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160 FEES EEE EE Sa SaSEEsuaes SETSSSEETSSEETTTEEET FTREEEETESEETTESEETS 


> SSEESRERRT ARERR ESA RRSE RESIS sesevatiass" 


Figure 7 


As I have already said, although Figure 2 shows the kind of model that is used in 
the analysis of mechanical vibrations, and Figures 3 and 7 look like the frequency 
response associated with such a model, you must not conclude that the procedures 
in this unit are useful only in the context of mechanical vibrations. On the 
contrary, they are widely used in such different fields as acoustics, meteorology, 
electrical network theory, and the design of measuring intruments; and this is by 
no means a complete list. In addition, Fourier series are mathematically useful in, 
for example, the solution of differential equations. 


So, in order to decide such different matters as the suitability of a given hi-fi 
amplifier for the reproduction of a particular tape recording or the relationship 
between the reading produced by a measuring instrument and the actual value of 
the measured quantity, we need the frequency response of the amplifier or the 
measuring instrument and a procedure for expressing the signal (i.e. the sound or 
other quantity of interest) as a sum of sinusoidal functions. The next section is an 
account of such a procedure. 


Summary of Section 1 
1. A function f is periodic if 
f(t) =f(t + nT) ier pt tor ee tt TS. 


where T is the period. 


See Unit 32, Partial differential 
equations. 
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2. The way that a mechanical system responds to sinusoidal inputs of different 
frequencies is called the frequency response of the system. 


3. A linear system is one that can be modelled by a linear diffierential equation. 
A linear second-order differential equation was defined in Unit 6 — this 
definition extends.naturally to higher-order equations. The superposition 
principle applies to any linear differential equation. 


4. If the frequency response of a linear system is known, the response to any 
periodic input, F, can be found provided that F can be expressed as a sum of 
sinusoidal terms, 


Pa? i454 + F,. 


The output, x, will then be the sum of the outputs, x;, to each sinusoidal input, 
F,, taken separately. 


5. The sum of a number of sinusoidal terms of different frequencies is a non- 
sinusoidal function. If the frequencies of the sinusoids are integral multiples of 
the lowest one, the fundamental angular frequency, w, say, then their sum will 
be periodic with period 27/q, i.e. the period corresponding to the fundamental 
angular frequency. 


2 Fourier series for periodic functions 


2.1 Introduction 


The aim of this section is to set up a systematic procedure for writing a given 
periodic function in terms of the sum of sinusoidal (or ‘harmonic’) components. 


From the last section you will remember that the sum of a number of terms like 
A, sinwt + A,sin2wt + A3sin3m@t + ---+ A,sinkaot (1) 


(where A,, A>, A3, ..., A, and @ are constants, with w > 0 and k a positive 
integer), will be a periodic function, f, of t, with period 27/@ which equals the 
smallest period of the term A, sin wt with the lowest, or fundamental, angular 
frequency. What we are aiming at is to start with a periodic function f and end up 
with a series like (1), with all the constants taking known values. Before we get 
down to this, I want to make three adjustments to the series (1), all of them aimed 
at making it more general, that is to say, capable of representing more functions. 


One adjustment arises from the recognition that what we have so far done with 
terms like sin wt and sin 2mt can also be done with terms like cos wt and cos 2ot. 
In other words, you can get a periodic function by adding up a number of terms 
like 

B, cos@t + B,cos2wt + B;cos3mt + -:: + B,cos kat. (2) 


To take account of this fact I shall represent the series which constitutes our initial 
function by the sum of a number of sine and cosine terms, like this: 


A, sinwt + A,sin2wt + A;sin3mt + -::+ A,sinkaot 
+ B, coswt + B,cos2wt + B;cos3mt + -:: + B,cos kat. 


This can be written more concisely as 


k 

)> (A, sinnowt + B, cos nat). 

n=1 

The second change I want to introduce removes another restriction which arises 
from the fact that the mean value of any function like A sinnwt or B cos nwt taken 
over any interval of length T is equal to zero. Any sum of such terms also has a 
zero mean value, like the function in Figure 1(a) of Section 1. Now the procedure 
I am about to show you is not restricted to functions with a zero mean. It can 
quite easily cope with other mean values. For example, the periodic function 
sketched in Figure. 1(b) of Section 1 has a positive mean value and our procedure 
is quite capable of dealing with a case like that. It can also cope with a negative 
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mean value. In order to allow for a non-zero mean, all we need to do is to shift 
the graph of our function up or down by the required amount, which simply 
requires an appropriate constant added to the series, so that our series becomes 


k 
f(t)}=M+ Y (A,sinnot + B,cos nat) (3) 


n=1 


where M is the new constant. 


The addition of the term M will not, of course, change the period of the right- 
hand side of Equation (3). As far as frequency response calculations are concerned 
(like those we did in Section 1), M is a cosine term with zero frequency and its 
contribution to the total response is worked out accordingly, as in the example 
below. 


Example 1 


Work out the effect on the response of the system described in Exercise 2 of 
Section 1 of adding a constant term to the input function f(t) so that the new 
input function becomes 


Sf (Hnew = 3 + sin 4t + 2sin 8t + 3 sin 12t. 
Solution 


The new constant will add an extra term to the response. It corresponds to a 
cosine term with zero frequency. From the graph in Figure 7 of Section 1 we see 
that the amplitude ratio is 1 and the phase shift is zero for w = 0. So the output 
term corresponding to 3 in the input is 3 x cos(0 + 0) = 3 and the total response 
will now be 


X(t)new ~ 3 + 1.2 sin (4t — 0.19) + 4sin (8t — 0.77) + 4.5sin (12t — 2.3). 


There is one more change I need to make to the series (3) so that it can represent 
a wider class of functions, and that is to extend it to an infinite series, i.e. 


f(t)}=M + Y (A, sinnat + B,cos not). (4) 
n=1 
The right-hand side of Equation (4) is known as the Fourier series for f(t). Now, 
we know that this series represents a periodic function with period 2z/m. But what 
we are trying to do is to find this series for any given periodic function. Can this 
be done? And, if it can, how do we find the appropriate values of the 
constants M, A, and B,? 


The answer to the first question, for practical purposes, is ‘Yes’. (We shall see what 
conditions the given function must satisfy later, in Subsection 2.3.) But for now we 
shall assume that our periodic function can be expressed as a Fourier series and 
get on with finding the constants. 


2.2 Finding the constants 


Let us see how to calculate the constants M, A,,; B,, Az, Bz, ... in Equation (4) 
from our knowledge of the function f(t) itself; this might be given, for example, in 
the form of a graph. We shall need to use the following integrals. 


Some useful integrals 
In all cases, m and n stand for non-negative integers and w = 2n/T. 


(5a) 


F/2 
sin nwt sin mot dt = | cos nwt cosmatdt = 0, n#m. (5b) 


~—Tis 
sinnwt cos mot dt = 0. (Sc) 


T/2 r 
sin? nwt dt = | cos’ not dt = —, (5d) 
—T/2 Z 


You can check these integrals 
using the Handbook. 
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Now to find M we integrate every term on both sides of Equation (4) with respect 
to t between the limits — T/2 and T/2. Hence 


T/2 Fi T7z T/2 
Jd = Mae + | A, sincotdt + | B, cos at dt 
= T)2 medt & UF - eat 8 F! = 1Ti2 
T/2 Ti2 
+ | A, sin 2mt dt + | B,cos2m@tdt+---. 
= Tie eee 


On the right-hand side, every term except the first one is zero according to (5a), so 
that 


T/2 


f(t)dt = TM 
T/2 


1 T/2 
therefore M =— f(t) dt. (6) 
a mae Se 


In other words, M is the mean value of S(t). 


We can find A, by multiplying every term on both sides of Equation (4) by sinawt 
as follows: 


f(t)sinwt = Msinwt + A, sin? wt + B, cost sinwt + A, sin 2otsin wt 
+ B, cos 2@t sinwt + °:° 


and now we integrate every one of these terms over t between the limits of — T/2 
to T/2: : 


T/2 T/2 T/2 
f(t)sin@tdt = M sinqwtdt + | A, sin? wt dt 
—T/2 —T/2 —T/2 
T/2 T/2 
+ | B, cos @tsinatdt + | A,sin2otsinot dt 
—T/2 sie 


Tiz 
+ | B,cos2wtsinwtdt+-:-. 
eet ie. 


If we compare all the terms on the right-hand side with (5a), (5b), (Sc) and (5d) we. 
see that they are all zero except the second one. Hence 


T/2 T/2 T 

f(@sineotdr = | A; sin’ wtdt = A, x — 

—T/2 —T/2 2 

SO 

f. T/2 

A, == f (t)sin ot dt. 

T J_rp 

Exercise 1 


By multiplying every term on both sides of Equation (4) by cos wt and then integrating 


each of the resulting terms over time from — 5 to x show that 


5 T/2 


B, == f(t)cos wt dt 
r Sal Bee 


where T is the period of f(t). 
[Solution on p. 23] 


All the other coefficients are found in the same way: by multiplying all the terms 
by the sinusoidal function which forms part of the same term as the required 
coefficient and then integrating over time from — T/2 to T/2. It follows, for 
example, that 


T/2 T/2 
f (t)sin 2mt dt = | A, sin? 2ot dt 
—T/2 —T/2 

so that 


7/2 


Z 
A,=— f (t) sin 2at dt. 
sere 
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Similarly 
Tf2 
B, =— f(t) cos 2et dt and so on. 
er 
Hence for any value of n 
2 at 
A,== f(t)sin not dt (7) 
PJ o7 
and 
> T/2 
B,==> f (t) cos not dt. (8) 
ig T/2 


By following the procedure in the box below, we can now find the Fourier series 
corresponding to some simple functions. 


Procedure 2.2 
To find the constants of the Fourier series 


f(t)=M+ y (A, sinna@t + B, cos nat) (4) 


n=1 


for a periodic function f(t) starting with a graph of f(t) against t. 


1. From the graph determine T, the smallest period of the 
function f(t), and hence the fundamental angular frequency w 
by putting w = 2z/T. 


T/2 
Write down the equation for f(t) dt in terms of t. 
=F /2 
1 T/2 
Find the constant term M = 7 f (t) dt. (6) 
T/2 


(This is the mean value of f(t) over one complete period.) 


Find the coefficients A,, Az, A3, ... of the sine terms by evaluating 


T/2 


2 
A, == f (t)sin not dt es 3) £7) 
T <2 f2 


(@ as determined in Step (1) above). 


Find the coefficients B,, B,, B3, ... of the cosine terms by 
evaluating 
| T/2 


B, =— 


f (t) cos not dt (n= 1, 2,3, ...) (8) 
T J-rp 


(@ as determined in Step (1) above). 


Example 2 


Derive the Fourier series for the saw-tooth function shown in Figure 1. 
f(t) 


| | 
| | | 
| | | 
| | 


| 
Figure 1. Saw-tooth function. Note that this function is not defined at t=..., —7, 7, 37,.... 
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Solution 

First, we find the period, which equals 27 because the graph repeats at intervals of 
2x. From this we find the fundamental angular frequency w to be 1 from the 
formula w = 27/T. It follows that the Fourier series has the form 


f(t)=M+ > (A,sinnt + B,cosnt). 
n=1 


Now M is the mean value and is given by 


Se eke h fy 
M =— (t) dt where — = 7. 
i ge Se ae 2 
For t between —7 and 7, f(t) ='t/x (except for the end-points, t = —z and t = 7, 


which do not affect the value of the integral). 
Therefore 


oe Lie 
= = SS ee = (), 
er | fe al 5 


(This was to be expected from the symmetry of the graph.) 
Again 
F/2 


2 : be ae 
A, == f (t) sin not dt = - | —sinnt dt 
T J _72 a pee 


1 Tt 
= =| t sin nt dt. 
Now, from the tables of standard integrals in the Handbook we find that 
t = 
fesin ata = — 7008 at + a2 sin at ee. 


In our case, a = n and for a definite integral we may omit the constant C, so that 


1 t os see 
A, = -3|-— — 00s 4-- sine 
Tl n n 


1 
— [—ncosnnx + 0 + (—2)cos(—nz) —0] 


TT 


ae 


_ —2ncosnn — —2cosnn 
ee 
Now cosnz = (—1)" so If you have not met this 
; before, note that 
a n n (=) = = 
A= (- ST. (== 
dese a and so on. 


(Remember that n is any positive integer.) 


So 
—] 2 —2 
A, =-, A, =—., Az; =—, A, =—., and so on. 
T Tt 37 4n 
Similarly 
2 ria 1 wr 
B, == f (t)cos nwt dt = — | t cos nt dt. 
T Jr No Jn 


Again, we find from the Handbook that 
ore 1 
[cos ade a 7 sin at + Pea at+C. 


Again, a = n and C is omitted for a definite integral, so that we have 


ii = ee 1 & 
B, = —| —sinnt + — cosnt = (). 
m™ | n n 


=e 
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Hence the required Fourier series consists only of sine terms since M = B, = 0; 


therefore 
comes: 
f®= > —(-1)""' sinnt 
n=1 7UN 
= ae sin 2t e sin3t sin4t = 
3 2 3 4 


Notice that the higher the frequency of a term, the lower its amplitude. 


Exercise 2 
How would the Fourier series derived in Example 2 above be affected if 


(i) the height of the ‘saw-tooth’ were doubled, i.e. the value of f(t) increased from —2 
to +2 as t increased from (2m + 1)x to (2m + 3)x (where m is any integer)? 
(11) the height of the saw-tooth remained the same as in Figure 1 but the period were 


doubled to 42? 
[Solution on p. 23] 


2.3 What kind of functions? 


What we have done so far is to show that if a function f(t) with period T can be 
expressed in the form 


f(t)=M+ Y (A,sinnot + B,cos nat) (4) 
n=1 
then the coefficients M, A,, B, can be calculated from the formulae (6), (7), (8). To 
make full use of this result we need to know under what conditions a given 
periodic function f(t) can be expressed in the above form. 


There is a theorem which gives a precise result of this kind; in order to state it we 
need two new mathematical concepts: continuous functions and piecewise 
continuous functions. 


A function is said to be continuous if its graph is an unbroken curve, that is if 
there are no gaps in it. For example the sine function is continuous, but the ‘saw- 
tooth’ function of Example 2, illustrated in Figure 1, is not, because of the gaps in 
the graph att = ..., —%, 3 3,.... 


A function is said to be piecewise continuous if its graph consists of continuous 

pieces; there must be only a finite number of pieces on any finite part of the 

domain, and the two pieces ending at each gap must have well-defined end-points 

even though these points need not themselves be parts of the graph. For example 

the saw-tooth function illustrated in Figure 1 is piecewise continuous; each piece 

has well-defined end-points even though these are not part of the graph. On the 

other hand the function 1/x is not piecewise continuous because the pieces do not The graph of 1/x is shown in 
have well-defined end-points at the gap where x = 0. the Handbook. 


Now we can state the theorem: 


Theorem 2.3 

If a periodic function f has period T and is piecewise continuous, 
and if its derivative is also piecewise continuous, then it can be 
represented by the Fourier series 


f(t)}=M+ Y (A, sinnot + B,cos not) 


n=1 
where wm = 27/T. The coefficients are given by 


T/2 


1 
M = rT ee (6) 


T {2 
A, == f(t) sin not dt (7) 
fe 
3 T/2 
B, == f (t) cos not dt. (8) 
bey 
The formula (4) is valid for all values of t at which there is no gap 
in the graph. 
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Exercise 3 
Derive the Fourier series for the ‘square wave’ function shown in F igure 2. 


| 
| | 
| | 
| | 
| | 
| | 


Sesser 


| | 
| | 
| | 
| | 
| | 
| | 
| | 
a 
| | 
| | 
| | 
| | 
| | 
| | 
| | 
| 


Figure 2. Square wave function. Note that the function is not defined at t = ae 
=1/2.172. 3/2 


[Solution on p. 23] 
Exercise 4 
Derive the Fourier series for the periodic function whose graph is shown in Figure 3. 


fo 


4 2 a 0 1 2 3 4 
Figure 3 
[Solution on p. 24] 


Exercise 5 


For the function shown in Figure 4, show that M = B.=0 and derive an expression for A,,. 


: f(t) : 

| | 

| | 

| | 

— ; | 

| | 

| | 

| | 

| | 

= — 0 h a t 

| | 

| | 

| | 

| = | 
ws 

| 

| 

| 


| 
| 
| 
Figure 4 
[Solution on p. 24] 


MST204 31.3 oe: 


Summary of Section 2 


1. The Fourier series for a function f(t) with period T is 


f(t)=M-+ Y (A,sinnat + B,cos nat), 
n=1 


where w = 2z/T. 


2. The constants in the Fourier series are given by See Procedure 2.2. 


1 T/2 
M =— t) dt 
T) 1” 


2. T/2 

A, == f(t) sin not dt (y =1, 2, 3,-.:.) 
T J_r 
a T/2 

B, == f (t) cos not dt i= 3 
T Jr 


3. (i) A function is said to be continuous if its graph is an unbroken curve. 

(ii) A function is said to be piecewise continuous if its graph consists of 
continuous pieces and if there are a finite number of such pieces on any finite part 
of the domain. The end-points of the pieces must be well-defined. 


4. Any periodic function which is piecewise continuous and whose derivative is Theorem 2.3. 
also piecewise continuous can be represented by a Fourier series. 


3 Some outstanding points 


3.1 Odd and even functions 


The object of this subsection is to show you how to use the symmetry of the graph 
of f (t) to save yourself work in calculating the Fourier series. 


In Section 2 we derived the Fourier series (or at least the values of the coefficients) 
for the four functions which are sketched in Figure 1. 


fi) 


| 


RAS MAUL ok Markie aaa | SMD ACALR Li 
S 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
I 


Figure 1. Periodic functions with symmetry about the origin. 


All these four functions have some symmetry about the origin. But the symmetry 
takes two different forms: in Figures 1(b) and 1(c) the graph is symmetical under 
reflection in the vertical axis; this kind of symmetry is described by the formula 
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f(t) =f(—t). This is the kind of symmetry which also applies to the function 

cos wt, and a function with this symmetry is called an even function. On the other 
hand, in Figures 1(a) and 1(d) the graph is symmetrical under a rotation about the 
origin and this symmetry can be described by the formula f(t) = —f(—t) which 
also applies to the function sin wt. A function with this symmetry is called an odd 
function. 


Now it is a fact that the Fourier series for an even function does not contain any 
sine terms (i.e. A, = 0); apart from the constant term it can only consist of cosine 
terms. Similarly an odd function does not contain any cosine terms (i.e. B, = 0) 
but only sine terms. 


This piece of knowledge can save you the trouble of trying to evaluate coefficients 
which you can immediately predict will be zero. Look back at the four Fourier 
series in Section 2 and check that they conform with what I have said about odd 
and even functions. 


Exercise 1 


Figure 2 shows the graph of a periodic function of t in the shape of a triangular waveform. 
It is drawn so as to be symmetrical about the line t = 0. It is an odd or an even function? 


[Solution on p. 25] 


If a functron of t has no symmetry about the origin then it is neither odd nor even 
and its Fourier series will contain both sine and cosine terms. 


Whether a function is odd or even or neither may depend merely on the choice of 


ee pepe Figure 2. Triangular 
origin. For example, in Figure 2, if the graph were moved to the right by a waveform. 


distance less than a quarter of a period, the graph would no longer have symmetry 
about the origin and the Fourier series would contain both sine and cosine terms. 


Exercise 2 


State which of the quantities M, A,, B, will be zero for the function which results when, in 
Figure 2, the graph is moved 


(i) one unit vertically downwards; 
(i1) one unit vertically downwards and 4 period to the left. 


[Solution on p. 25 | 


Exercise 3 
Derive the Fourier series for the function whose graph is shown in Figure 2. 


(You can take advantage of the symmetry when doing the integrations for the terms which 
are not zero by noting that, for example, 


[fla 2 food 


because f(t) is even.) 
[Solution on p. 25] © 


3.2 The use of Fourier series as approximations 


As we have seen, a Fourier series generally has an infinite number of terms. 
Clearly, in any numerical calculation we can cope with only a finite number of 
terms and the question therefore arises: how many terms do we need to get a 
good approximation to the original function? Of course, it all depends on what we 
mean by a ‘good’ approximation; and this must depend on what we want the 
series for. No general specification can be given, but some ‘feel’ can be obtained by 
looking at some functions like those we dealt with in Section 2 and comparing 
their graphs with those obtained by adding up an increasing number of terms of 
the Fourier series. You may be surprised at how close even quite a small number 
of terms will take us. 


By way of an example, we can start with the saw-tooth waveform, the first Fourier 
series we derived. Figure 3 shows the results of using from one to six terms of the 
series. 
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fd) fi) 


using one term 


using two terms 


f() f(t) 


using three terms using four terms 


f(t) fp) 


using five terms 


using six terms 


Figure 3. Graphs of partial sums of the terms of the Fourier series for the saw-tooth 
function. 


Figure 4 shows a triangular wave with, superimposed, the fundamental component 
(curve I) and the sum of the first two components (curve II) of the Fourier series. 
Neither curve is very far removed from the original function. 


Figure 4. Triangular wave with the first two partial sums of the Fourier series 


8 1 1 
a= 73 (cos t + gz 008 5 + ry tha St + °°°). 
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Clearly, then, even if a Fourier series is ‘truncated’ (i.e. cut off) after a finite 
number of terms it can still give quite a reasonable approximation to the original 
function. 


Summary of Section 3 


I. A function which exhibits symmetry about the origin such that f(t) = f(—t) is 
known as an even function and its Fourier series contains no sine terms (A, = 0). 


2. If the symmetry about the origin is such that f(t) = —f(—t) the function is an 
odd function and its Fourier series contains no cosine terms (B, = 0). 


3. Fourier series are used as approximations to the original functions by using a 
limited number of terms. For many purposes good approximations can often be 
obtained with quite a modest number of terms. 


4 Extending the scope 


4.1 Functions of variables other than time 


So far the main applications we have had in mind have been those concerned with 
variations of time. It is quite possible, by analogy, to use Fourier series in cases 
where the independent variable represents a quantity other than time (e.g. space). 


Suppose we have a function of some other, unspecified variable, x. Suppose further 
that the function repeats itself at intervals of length L along the x-axis, as shown 
in Figure 1. We can derive the Fourier series exactly as we did before, provided 
that we substitute x for t and L for T. We shall then finish up with a Fourier 
series in x, as follows. 


Fourier series in x 


f(x)=M+ > (A,sinnwx + B,cos nox) 


n=1 


4.2 Some non-periodic functions 


So far, we have dealt only with periodic functions. It is possible, however, to 
extend the techniques of Fourier series to some non-periodic functions. This aspect 
of the subject is important mainly as a mathematical technique and this is how it 
will be used in the next unit. It is of very limited use for solving the kind of 
‘frequency response’ problem which we discussed in Section 1 because that kind of 
problem assumes a periodic input. The fact that I refer to ‘some’ non-periodic 
functions suggests that there is a restriction on the type of function we can cope 
with. The restriction is that the function must have a finite domain: it is defined 
over only a finite part of the t- or the x-axis. 


fix) 


L———+ 


Figure 1. Periodic function of x. 
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The procedure for finding a Fourier series to represent such a function really 
amounts to extending the definition of the function to the whole of the t- or x-axis 
in such a way that the extended function is periodic. Then we can find the Fourier 
series for the extended function but use it only over the range for which the given 
function is defined. Let us take a simple example, the function F(x) shown in 
Figure 2. 


Here F(x) = 2 — x and the domain consists of the interval of the x-axis between 0 
and 2. The trick is to incorporate this function in a periodic function for which the 
Fourier series can easily be found. One way of doing this is shown in F igure 3. 


I have drawn dashed lines to define an even function f(x) which is periodic with a 
period 4. Having got the Fourier series, we must remember to use it only over the 
range, 0 < x < 2, of the actual non-periodic function we started with. The Fourier 
series is like the one we derived in Exercise 3 of Subsection 3.1. Substituting x for 
fit is 


Figure 2 


1 3 1 e 
ea oes a 9 (for all x). 


g 
a 
LOSE et BO 


This can be used to represent the function F(x) shown in Figure 2 provided that 
we restrict ourselves to the domain 0 < x < 2, ie. Figure 3 


F(x) = 2 ea he BU et ean ) Qe x22 
x)=2—-X= —z (cos — + ~cos—— + ~cos—— + ::: 
2a ee co 
Exercise 1 
Using the first four terms of the series derived above, investigate how accurately this partial 
sum represents 
P(x) = 2 —x forO<x <2 


by substituting x =0, x=4, x=1, x=2. 
[Solution on p. 25] 


Exercise 2 


Again using only the first four terms, verify that the formula for F(x) given above cannot be 
used to represent the function 2 — x when x = 3. 


[Solution on p. 25 | 


Exercise 3 


Sketch the graph of a periodic function which could be used to represent F(x) as defined in 
Figure 2 by means of a series containing only sine terms. (Do not calculate the coefficients.) 


[Solution on p. 25] — 
f(t) 


Exercise 4 1- 


Derive a Fourier cosine series with M = 4 and period 4 to represent the non-periodic 
function shown in Figure 4, for the range 0 <t < 1. 


[Solution on p. 26] 
0 

Summary of Section 4 Figure 4 

1. The procedure for determining Fourier series for functions of time described in 


earlier sections can be used for functions of any variable x, with a period L, 
subject to substituting x for t and L for T. 


2. A non-periodic function whose domain is a finite interval can be represented 
by a Fourier series by deriving the series for a periodic function of which the given 
function forms a part and then restricting the independent variable to the range of 
values of the variable for which the non-periodic function is defined. 
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5 End of unit problems 


Please refer to the study guide for this unit before you begin work on this section. 


Problem 1 
Derive the Fourier series for the periodic function one complete cycle of which is shown in 
Figure 1. 


Figure | 
[Solution on p. 26] 


Problem 2 
Figure 2 represents the frequency response of an idealized device. 


output amplitude 
input amplitude 


Figure 2 


Show that, in response to any periodic input, the output of the device will be an exact 
scaled version of the input, shifted by an amount 1 along the time axis. 


[Solution on p. 26 | 


Problem 3 
(i) Derive a Fourier sine series of period 2/ to represent the non-periodic function whose 
graph is shown in Figure 3 over the range 0 < x </. 


(ii) Assuming that we may neglect terms in the series having an amplitude equal to or less 
than 1° of the amplitude of the first term, how many terms of the Fourier series are 
required to represent f(x) to this degree of approximation? 


[Solution on p. 27 | 


Figure 3 
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Appendix 1: Solutions to the exercises 


Solutions to the exercises in Section 1 


1. (i) The fundamental angular frequency is 2 radians per 
second. 

a pg et 
(ii) The period is ae seconds. 

2. (i) Since the angular frequencies of all the components 


are integral multiples of the lowest one, i.e. w = 4, the period 
T, of the function f(t) will be 27/4 = 2/2 seconds. 


(ii) The first component has angular frequency w = 4. At 
this frequency the relative response, or amplitude ratio, is 
found from Figure 7 to be about 1.2, and the phase shift 
about —11°, or —0.19 radians. Hence the output 
corresponding to the term sin 4t is about 1.2 sin (4t — 0.19). 


9 


Similarly at @ = 8, the amplitude ratio is about 2 and the 
phase shift about —44° or —0.77 radians. The output 
component corresponding to the term 2 sin 8t is about 

2 x 2sin (8t — 0.77) = 4sin (8t — 0.77). 


Again, at m = 12, the amplitude ratio is about 1.5 and the 
phase shift is about — 132° or —2.3 radians. The component 
of the output corresponding to the term 3 sin 12t is about 

3 x L.5sin (12t — 2.3), i.e. 4.5 sin (12t — 2.3). 


Hence, by the superposition principle, the total response of 
the given linear system to the input f(t) is 


x(t) ~ 1.2 sin (4t — 0.19) + 4sin (8t — 0.77) 
+ 4.5sin (12t — 2.3). 


(Your answers may differ slightly because the graph is hard 
to read accurately.) 


Solutions to the exercises in Section 2 


T/2 if2 


1. f(t)cos wt dt = 
~T/2 ~T/2 


M cos wt dt 


fips 
+ | A, sin wt cost dt 
=a Bl? 


T/2 
+ | B, cos? wt dt + -°: 
sed (Pe 


All the terms on the right-hand side except that containing 
B, are equal to zero. Hence 


Lis T/2 T 
| f(t)cos at dt = { B, cos? wt dt = B, x — 
=F —T/2 2 
Therefore 


a T/2 


B, == f (t) cos cot dt. 
T Jen. . 


2. (i) The period would still be 2z, and so w = 1 as 
before. Also, from the symmetry of the graph, M = 0 as 
before. So the new Fourier series has the form 


f(t)= > (A, sinnt + B, cos nt). 
n=1 
Now, for —1 <t < oe = 2t/m, so 
T/2 
=— =| = — ate nt dt = 


T/2 7 


St gee 
=| *t sin nt dat. 
a a 


This is just twice the original value’so 


4 
A. a —(-1)"*), 

mm 

Also 
T/2 
— == [ ane nt dt 
T/2 7% 
and again, since this is just twice the original value, B, = 0. 
So the new Fourier series is just twice the previous one: 
sin3t sin 4t 
3 4 


(ii) The period is doubled, so T = 4x and w = 4. The 
Fourier series now has the form 


sin 2t 


fie) = 9 sine ~ 


nt 
ft)=M+ y la, »sin + B, cos" 


n=1 


For —2n<t< 2n, f(t) =.t/2%, SO 


| 2n t 1 t? 22 
4n J_>, 21 ol Sage See 


(as expected, since the mean value is still zero). 


SS coat ee aes nt 
AS | =-- ee ie t sin — dt 
On "sin! - ae oe 


1 = : nt - 4 . nt 
= —.~| —-tcos— + —sin— 
4n7} on pe sie a Se 


2 
= [—2ncosnn + 0 — 2ncos(—nz) — 0] 


4n?n 
2 
= —(—1)"*! 
7Nn 
4 2n 
= t cos — dt 
i An? ie a 


So the new series is 


fo = 2 (sina) Sat 5 S02)... 


and the only difference is the changed value of w. 


3. From the graph, the period is 2. So w = 2n/T = 1. 


Here, 
Liz de) 


Be [fede i dt 


=< T72 


since over the rest of the cycle f(t) = 0. This time the mean 
value is clearly not equal to zero since no part of the wave- 
form is below the axis. So 


ee ee — 1/1 ; 1 
M =— t)dt = a =-|~+-—]=— 
ee A ee a= 3; 7 ee: 
T/2 1/2 
A, == f(t) sin not dt = 1 a 1 sin nat dt 
T Je 4p 


1/2 


1 1/2 
—-— COs nt 
nn —1/2 


| 
| - 
aa 
Q 
ro) 
n 
SS 
| 
(@) 
o) 
n 
—————. 
bil 
a 
[_—— 
SEA OS 
I 
ee 


2 rti2 1/2 
B, == { f(t)cos nat dt = [ cos nt dt 
TJ 72 wo ade 
E _ “| s =) 
= | —sinnazt = —| sin— — sin 
_ [nn ~1/2. NT 2 2 
y) : nt 
= — sin—. 
nt 2 
2 
Thus f= 35 = — 
2 us 
2 
B, =~ sinn = 0 
B 2 ian 2 
— § ed 
weed ee 3x 


and so on. 


Since the values of B,, are the coefficients of the cosine terms 
the required Fourier series is: 


- 


cos3mzt cosSzt  cos7zt 


5 7 

The brackets contain an infinite series of cosine terms (the 
sine terms all being zero), and again, the higher the 
frequency the lower the amplitude. 


2 
4. From Figure 3, T= 4= au 
@ 


LZ 
f(t) ==+-—|coszt — 
a 


252K 
H = = =, 
ence @ 4 5 


: 1 
The mean_value = M = = 


9 2 
A == { t? sin noot dt. 
2 
The integral tables in the Handbook give 
| an i? Dike. 
t“ sinatdt = | —-— + —] cos at + —sinat + C. 
a 2 a 


Here, a = nw so 


2 
ee aos | pareniote ae 
| 2 :} cos ne) | 
ee [2 sin 2nw@ — (—2si 2nw 
sin 2n — - 
jar sin (—2nq)) ] 


= 0. 
+2 


| 1 
Also B, = 3 | t? cos not dt. 


—2 


From the Handbook 


- Dear 2 Be 2t 
t“ cos atdt = |— — —] sinat + —cosat+C. 
eee a 


MST204 31 Solutions 


Now a = na, as before, hence 


tLe 2 re A = 
2s a sin not S79 fat = 
1 Z : 
ae (4 — os sin 2n@ 
2 : 
_ [4 - za] sin(—2n0) 
n 


; | 
ee [2 cos 2nw — (—2 cos (—2nw)) ]. 
no 


But w = 7/2, so that sin2nw = sinnnz = 0, and 


Fypceeatas al Pe 
ees cosnm = ——-(—1)". 
oe nn UR 


6 16 
B,=— x (-l= -s 
1 n2 ( ) m2 
4 4 
ee ee 
16 16 
i = oni (—1)? = oa and so on. 
So 
f(t) a ost a t + —cos urea 
= -=-—-— |cos~t ——cos at + ~cos— 
a oe 2 4 
4 162 (-1)" ono 
Soe ook 
eS 2 =z 0085 


5. From Figure 4 the period, T = 6h. 
So @ = 22/6h = 2:/3h. 


T/2 


1 3h 
M = 7 ae = a} flea 


1 =—h h 3h 

= — — hat + | tdt + | hat 
al —h h 
1 

= Gy (—2h* +0 + 2h*) = 0. 


(In fact we can see directly from the graph that the mean 
value is zero.) 


T/2 1 3h 
A, == f(t) sin not dt = — f(t) sin not dt 
io T/2 3h J 3, 
1 =i h 
all = hsinnestde + | tsin not dt 
3h \ J — 3p —h 


a3h 
_ | hsin nwt it 
h 


Using the integral tables in the Handbook to find the second 
integral 


1 t 
= sy) ~ pcos" =e 


1 1 3h 
= —— COs nat 
Al nw | 


1 
eo [cos (—nwh) — cos (—3nw@h) — cos 3nwh + cosnah | 
no 


1 h 
sin not 
no? is 


1 
>| — hcosnwh + —sin nah — hcos (—nawh) 
3na@h nw 


1 
— —sin (— na) | 
nw 
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The terms in cos nwh cancel out, leaving 


2 1 
ss ee nah — hcos Siebel. 
~ 3ncoh 
But @ = 2/3h, so 
2 i 3h..nt |. 
A, = —|—sin— — hcosnn 
nt \nt 3 
T/2 1 3h 
Bios 7 = not dt = 3h wl cos not dt 
1 


—h h 
-a( J — hcosnat dt + t cos not dat 
—3h eee | 

3h 


aa hcos nat i : 
h 


Using the integral in the Handbook for the second integral, 


| ee I 
B, = =| ——sin not 
3 n@ Soe 


fe —— sin nest + 
<j Si 
3h 


1 1 3h 
+ —| —sin nwt 
3 i | 


1 
=. [— sin (—noh) + sin (—3nq@h) + sin 3n@h —sin nh | 
nw 


1 1 
+ ——| hsinnwh + —cos nwh + hsin (—nah) 
3nah n@ 


1 
=——cps (— no) 
n@ 


a0 F 
Hence M = B, = 0 and 
2h { 3 nt ; 
= — |—sin— — cosnnz 
~ nan\nn 3 


Solutions to the exercises in Section 3 


1. The function is even since f(—t) = f(t), and its Fourier 
series would consist of a constant term (since its mean value 
is non-zero) and cosine terms. 


2 (i) fa 


This will be an even function with zero mean value i.e. 
M= AL =A. B,, # 9. 
(il) fo 
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This will be an odd function with zero mean value i.e. 
M= 8, = 0, A, #0. 
3. From the graph, T = 4, so w = 2n/4 = 1/2. 


|e ig ee ee 
M=;| sod=F[ soar 


Between 0 and 2, f(t) = 2 —t, so 
Mee sf Go yas : of oy Sy 
pS ee oe 
Since the given function is even, A, = 0. 


2x2 
Pace 


2, 
| (2 2 feng 
0 2 


- t t 
= I [20s _ reos "| dt 
= s yee Oe nit 4 nit [ 


2 nt eae: S nn? 2 


Tsar Geer 
nim 


4]. : 1 1 
= —/sinnz — sinna — —cosnz — sin0 + 0 + —cos0 
nt nit nt 


4 
= 522 | — COS nT). 


8 8 
Thus B} =—, B,=0, B3= O72? and so on.. 
7 
Hence the required Fourier series is 


fie i+e ee a. 2 ee 
Bd wee, ee a ee 


p 


Solutions to the exercises in Section 4 


1 
bh 2 Gj)? sitelts + —] ~ 1.933. 


Geo) 
F(S) = 15 and 1+ = (0.70711 — 0.07857 — 0.02828) 


~ 1.487. 
F(1) =1 and 1+ — 5 0+0+0)=1 
F(2) = 0 and 1+ ~ 0.067. 
= -5 25 


So that even with only a few terms the series gives an 
approximation to F(x) = 2 — x with a maximum absolute 
error of 0.067 for the stated values of x. 


8 
2. 2-3=-1 and 14-040 +40)=1. 
m4 


This confirms that outside the stated domain the series does 
not represent the original function. 


3. The required graph is shown below. 


f(x) 


Note that the full line represents the function G(x) = — 
0 < x < 2, so that 2 would need to be added to the sine 


3 
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series representing this graph in order to obtain the 
required representation of the function F(x) = 2 — x, 
jee Sas 


4. We need to incorporate the function shown in Figure 4 
in a periodic function whose Fourier series will then 
represent the given function over its range of values of t, i.e. 
C= f< t. 


Exercise 3 in Section 2 provides a suitable function which, 
with a little adaptation, will give us what we want. We first 
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draw the appropriate diagram in which only the given 
function appears in full lines. 


This is like Figure 2 in Section 2, except that the period is 
different. Here T = 4, so that w = 2/2 whereas in Section 2 
T = 2 and wm = 72. This, however, is the only difference so 
that we can use the same series as we obtained then, 
changing only the frequencies of the terms. We therefore put 


Lt ee 
f(pj= 5+ is cost a 


. ae eg lor 02 4-= I 
— COS — lf — —Ccos— eee : 
5 008-5 ee or 


You can check by substitution that this series works over the 
stated range 0 < t < 1. With a ‘jump’ discontinuity at t = 1, 
the series will give the average value of f(t) for the vertical 
portion, i.e. for t = 1 the series gives f(t) = 5. 


Appendix 2: Solutions to the end of unit problems 


1. From the graph in Figure 1, T = 2l, so m = 22/21 = m/l. 
The graph is symmetric about the t-axis so M = 0. It is an 
odd function so the Fourier series has sine terms only, i.e. 
B, = 0. 


2 T/2 1 l ; t 
A, = — f(t) sin not dt = “| f(t) ea 
T Jr a ee 
1 ay aE 
== 3 2 t 
So { Shain a + | hsin™ at 
l ~E- l aoe l 


z l l 

h([{ 1 rt oe ae ie l nit |2*4 
= —<{/ — cos — + | ——cos— 

l (| na l -he nt l rd 

oh nt ae [nt l 
= *| cos" (-5 + 4]] — cos (> al] 
nr | 1 nt | | 
~ 00s ("= (5+ d}) + cos ("* (5 — a}} | 
d 


2h ; "ap a 
=—| cosnn ie a 
nt ae 


1 

~ —-—] — cosnn 

l 

Using the trigonometric identity from the Handbook 


2 


sina sinB =4cos(a — B) —4cos(« + 8), 
4h . nn . nad 


A, = — sin—sin —. 
ene - l 


Hence 


© 4h d t 
ti 22 _ sin = sin = sin — 


2. Let the input be f(t) where 


f(t)= M+ A,sinowt + A, sin 2ot + A; sin 3ot + °°: 
+ B,coswt + B,cos2mt + B;cos3mt+-::-. 
Now the amplitude of each of these terms will be multiplied 
by py in passing through the device since the same multiplier 


applies to all frequencies. (For this purpose M counts as a 
cosine term with zero frequency.) 


Each term will also experience a phase shift given by 
b= —TO. 

For M, dy =tx0=), 

for A, sinwt, @; = —Ta, 


for A, sin 2ot, od, = —2TH and so on, 
and similarly for the cosine terms. 
Hence, if the output is x(t), then 
x(t)= uM + pA, sin(@mt — wt) + WA: sin (2at — 2m@t) + °°: 
+ pB, cos(wt — wt) + wB, cos(2mt — 2wt)+ °°: 
= u(M + A, sinw(t — t) + Az sin2@(t —t)+ °°: 
+ B, cosw(t — Tt) + Bz cos 2m(t — t) + -°*). 
This is an exact scaled version of the input f(t) where (t — Tt) 
is substituted for t, i.e. each term is shifted by an amount t 
along the t-axis. The ideal frequency response shown in 
Figure 2 of Section 5 is not, in practice, attainable. You may 


be interested to see how near to the ideal a practical 
measuring device in common use gets. 


amplitude ratio (%) 


BO UF Oe ee 


frequency ratio = applied frequency 
natural frequency 


frequency ratio 


radians 
o( ) 02° 0.4: $:6-- 02 et Ss 


0 


— 1/4 


—7/2 


MST204 31 Solutions 


The above figure shows the frequency response graphs for a 
galvanometer on a commercially available voltage recording 
instrument. You can see that up to a frequency of about 0.7 
times the natural frequency of the galvanometer, the actual 
relative response is within +5% of the ideal, and the phase 
shift graph is quite close, too. For frequencies higher than 
about (0.7 x natural frequency) the relative response departs 
very markedly from the ideal. Hence, in order to get a true 
recording of a periodic input from this instrument, the terms 
of the Fourier series of the input with frequencies higher than 
(0.7 x natural frequency of galvanometer) must be negligible 
or zero. 


3. (i) The question asks for a sine series, i.e. a series 
consisting only of sine terms; this means that we must find 
an odd function of which the given figure is a part. Such a 
function is shown below. 


Here the ‘fictitious’ parts of the periodic function are drawn 
in dotted lines. For this periodic function, M = B, = 0, and 
L= 21,30 Gat @ = 22/0 = #1. 


9 l 
A, == ah f (x) sin nwx dx 


1/2 ee See 
zo a (- x — 1)sinnawx dx + “| — x sin n@x dx 


—1/2 


Pes fe —x + l)sinn@x dx 
l Jy l 


22a E 
= —; cosn@x — 
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4/2 


ae 

7 Sin nw@x + —cOSnwWx 
2m now a 
1/2 


sin nox 


(60) 
bi 1 
160) n-@ —1/2 


n 

x 1 l : 
—cosnax — =—sinnwx — 7 008 nwx 
nw n*w? 


“fol a oo 


1 
+ lcos (—nal) + —sin(—nol) — I cos (—nal) 
nw 


9g 
2d 
os 
2d 
+P 
2d 


l ee nol : 1 Ss nol ie nol 
es Se ees 1 ee ee ae 
jain, ae ae Mee, 


not a. eae 
a + lcosnwml ——sin nal — lcos nal 
2 no . 


l + now See nol 
yo 2 


oe ae 3 ay 
~ Pho 2 ee ee 2 


+ (l1—1+1-—I)cosnal 


| 1 1 1 1 | nal 
ee 
nw no no no 


+ | ——— — —]sin nol 
nw@ nao 


2d 


I7n?@? 


(4 sin — 2sin no} 


Zz? 
1 A 
But @ = T so sinnal = 0 and 
8d . nt 
A, Bas a 
8d 8d 8d 
Thus A, = -—5, A, = 0, 4, = ——~, A : ler 
us A, ar Aa 3 9,2” 44 0, A; 7572 and 
SO ON. 
Hence the Fourier sine series is 
1X ee ees 3 
ify= sin — ee a et pet, 
for0<x<l 
or 
8d (oe) 1)"*1 3 7 1 
TiS a Lea Oo ain for 0 ee 8S lL. 


ee . = 8 
(ii) The amplitude of the first term is ok 
Tl 


The ratio of the amplitude of the nth term to that of the first 


1 
(eri is = 
(2n — 1)? 


: 1 
For this to be greater than 1% we want (On — 1? > 100 1.€. 
(2n — 1)? < 100 or 2n — 1 < 10. This means for n < 5 we 
must include the terms in the Fourier series, so five terms are 
required. 
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Appendix 3: Television programme notes 


You should read these notes before watching the television programme ‘Fourier 
Analysis’. 


A.1 Introduction 


The television programme for this unit is about measurement and recording. It 
describes a procedure for predicting the response of a recording device to any 
periodic input that can be represented by a Fourier series. This procedure can be 
used for instruments such as audio amplifiers, record players or 
electrocardiographs. In the programme we use it to test a galvanometer which 
forms part of an ultra violet light recorder. 


In order to use the procedure we need two pieces of information: 
(1) the frequency response of the recording device; 
(11) the Fourier series of the periodic signal to be measured or recorded. 


Since the Fourier series has an infinite number of terms it must be truncated to be 
used in practice, but enough terms must be left to constitute a good 
approximation to the original waveform. As I said in Subsection 3.2, there is no 
hard and fast definition of a ‘good’ approximation; the programme models the 
triangular wave, which serves as the main example, as the sum of only the first 
three terms of the relevant Fourier series. (You can get an idea of how good this 
approximation is by looking at Figure 4 in Subsection 3.2 which shows the sum of 
the first two terms in this Fourier series compared with the complete waveform. ) 


The programme also shows how the relative response of the galvanometer is found, 
by providing a constant-amplitude sinusoidal input over a range of frequencies and 
noting the variation in the output amplitude. 


A.2 A note on frequency 


In the programme, frequencies are quoted in cycles per second because 
instruments are calibrated in this way. However, in this unit we have worked 
exclusively in angular frequencies (radians per second). The connection between 
these quantities is quite simple. 


If jf = frequency; in cycles per second (or Hz) 
and w =angular frequency, in radians per second, 
then 


and 
= 27). 


Example 1 


(i) The fundamental frequency of a periodic waveform is 30 cycles per second. 
Calculate the corresponding value of the angular frequency. 


(11) The component with highest frequency in a truncated Fourier series has an 
angular frequency of 600 rads" *. What is its frequency in Hz? 
Solution 
(i) @ =2nf~ 188.5 (in radians per second). 
600 


Gi) f=— = 5 4 
Ps Cees 


The relative response, i.e. 
output amplitude 
input amplitude ’ 

is called the amplitude 

response in the television 

programme. 


The SI unit of frequency is the 
hertz: 
Lie = bees = 


f and w are related to the 
period, T (in seconds), by 
1 2n 


ot 
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_A.3_ Relative response and phase shift 


As you saw in Subsection 1.2, the frequency response of an instrument consists of 
its relative response (i.e. output amplitude/input amplitude) and its phase shift. 


The relative response of the galvanometer in the ultra violet light recorder is 
shown in Figure 1(a). For good reproduction, the graph of the relative response 
should be as nearly as possible parallel to the horizontal axis over a range of 
frequencies. This range should include all the frequencies in the truncated Fourier 
series representing the signal to be recorded. In this case you can see that the 
graph is flat up to a frequency of about 90 cycles per second. 


output amplitude 
input amplitude 


100 ‘frequency 
(a) (cycles per second) 


frequency 


phase shift 
(cycles per second) 


100 


0) 


(b) 


Figure 1. The frequency response of the galvanometer tested in the television programme. 


The phase shift arises because the recording device takes time to respond to an 
input signal so that the output is not simultaneous with the input; each 
component of the Fourier series representing the input is delayed by an amount 
which depends on its frequency. Even if all the output amplitudes are correct, they 
will not add up to a good reproduction of the input if the instrument introduces a 
time shift of one component relative to another. In other words, the time shift 
should be the same for all components. To see what this means in terms of the 
frequency response, consider just one sinusoidal component, f; = F; sin wt, say, of 
the input, f(t). Then the corresponding component of the output, x(t), will be 

x; = X;sin w(t — tT), where t is the delay introduced by the instrument. But in 
terms of the phase shift, ¢, the output component is x; = X;sin (wt + ¢) so 


o = —Qt. (1) 
We want the time shift, t, to be constant for all values of w, so Equation (1) shows 
that the graph of ¢@ against w should be a straight line. The relevant graph for the 
galvanometer is shown in Figure 1(b), and you can see that the frequency range 
over which this graph approximates a straight line is greater than the frequency 
range for which the relative response graph is flat. So provided we stay within the 
frequency limits imposed by the relative response, we shall automatically satisfy 
the phase shift requirements as well. For this reason the programme concentrates 
on the relative response. 
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Example 2 


Assume that the triangular waveform shown in Figure 2 can be represented by the 
sum of the first three terms of its Fourier series; estimate the highest fundamental 
frequency at which such a wave will be reproduced by the galvanometer used in 
the programme. 


Solution 


From the solution of Exercise 2(i) in Section 3, M = A, = 0. Modifying the 
solution of Exercise 3 in Section 3, we find that the series is 


fith= = (cos at + = c0s 3at + 55008 Sat) 
where w = 27/T is the fundamental angular frequency and K is the height of the 
triangular wave. 
Now, in order to stay within the flat region of the relative response graph (and 
* hence also the linear region of the phase shift graph), we must have 
So < 22 x 90 


2n x 90 
oS 

5 
The fundamental frequency must be less than about 113 radians per second. 


~ 113. 


Now watch the television programme ‘Fourier Analysis’. 
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